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The scale of low energy cc and bb cross sections on light hadrons is of great importance to searches 
for the quark gluon plasma using the heavy-quarkonium suppression signature. Unfortunately, little 
is known about these near-threshold cross sections at present, and recent theoretical estimates span 
many orders of magnitude. Here we use experimental data on the four observed closed-flavor heavy 
quarkonium hadronic decays tjj' — > nnj/tp, T' — > nnT, T" — > nnT and T" — > mrT' , combined 
with simple models of the transition amplitudes, to estimate the pion scattering cross sections of 
cc and bb mesons near threshold. Specifically we consider the closed-flavor reactions irj/ip — > Trip', 
7rT — > 7rT', 7tT — > 7rT" and nT' — > ttT" and their time-reversed analogues. Our results may be 
useful in constraining theoretical models of the strong interactions of heavy quarkonia, and can be 
systematically improved through future detailed studies of dipion decays, notably ip ' —* irnj/ip and 

Til 

PACS numbers: 12.38.Mh, 13.25.Gv, 13.75.Lb, 25.75.Dw 



INTRODUCTION 

One signature proposed for the identification of a 
quark-gluon plasma (QGP) is a suppressed rate of for- 
mation of charmonium bound states, due to the screen- 
ing effect of the plasma on the linear confining potential 
that normally binds a cc pair [lj. In the presence of 
this screening, cc pairs produced within the plasma pre- 
sumably separate into open charm mesons, so that the 
formation of a QGP would be signaled by a decrease in 
the charmonium production cross section. 

This attractively simple picture becomes more compli- 
cated if the dissociation cross sections of charmonia on 
light hadrons are not small. In this case the charmonia 
produced during a heavy-ion collision may rescatter into 
open-charm final states as the result of interactions with 
the "comoving" light hadrons produced in the collision. 

Due to the importance of these cross sections for QGP 
searches, many calculations of near-threshold scattering 
cross sections of light hadrons on charmonia have re- 
cently been reported. The methods applied include a 
high energy color-dipole scattering model 13, quark mod- 
els HUE!, meson exchan ge m odels 0, 013 > an d most 
recently QCD sum rules |^llCJ|. 

Under certain simplifying assumptions one may relate 
experimentally known hadron decays to heavy quarko- 



nium scattering cross sections. Here we consider the dip- 
ion decays of heavy quarkonia below open-flavor thresh- 
olds, which have been observed in four cases, tp' — > 
7T7rJ/> EL El EH, T' ttttT [HEIEIEI T" -> ttttT 




and T" -» ttttT' [HEl El El El 



These transitions are quite weak, with partial widths of 
only ca. 100 keV for cc and 1—10 keV for bb. They have 
been observed only because the total widths of the ini- 
tial heavy quarkonia, which lie below their open-flavor 
thresholds, are also very small. 

We can use these three-body partial widths to esti- 
mate pion scattering cross sections of the corresponding 
heavy quarkonia, since these processes are described by 
the same invariant amplitudes. The pion scattering cross 
sections are given by 
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the limits are ti = —2[Et Ti E. K2 TPiPf — m V\i an d frame- 
dependent quantities such as E 7Tl and pi = \pi\ are un- 
derstood to be evaluated in the cm. frame. 

The corresponding dipion three-body partial width is 
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FIG. 1: Dalitz plot of kinematically allowed regions for the estimated from J/v , assuming constant amplitudes, 

decay ip 1 — > n + -k~ J / tp and the reactions ir + J/ip — > 7r + i/> ' and 



where s and t (after crossing) become s = m\^ M and 
t = m 2 ni7r2 . Our Eqs.(l,2) correspond to Eqs. (37.30,21) 
of the 2002 PDG [Sj. 

Our convention is that Alb is higher in mass than M a , 
so the dipion decay Mb — * 7T7rA4 a is energetically allowed 
but M a — > 7T7rAlb is not. The initial and final mesons 
in the cross section formula Eq.(l) may be (M a ,Mb) as 
shown, or they may be transposed to (Mb, Ma), since the 
same invariant amplitude A is involved. As an example, 
the decay if>' — > Tr + ir~J/ip is related to both ir + J/ip — > 
7r + V>' and ir + ip' — > ir + J/tp. 

Since the squared invariant amplitude ( |-4| 2 ) is sam- 
pled in different kinematic regions by the decay and the 
cross sections (see Fig.l), assumptions regarding the form 
of ( |.4 | 2 ) are required to relate these various processes. 
In the following we shall obtain results for the cross sec- 
tions given three simple models of the invariant ampli- 
tudes. 

CONSTANT A APPROXIMATION 

As a first approximation we neglect any dependence of 
( | 2 ) on kinematics, and simply treat it as a constant. 
In this case Eqs. (1,2) imply a simple relation between the 
pion cross section and the dipion partial decay width, 

_ 167r 2 Af b 3 Pf _ p f 

OTrMa^-rrMb — 1 Mb — >TrirM a 'a b — Co * 

A D Pi Pi 

(3) 

where Ajj is the area of the Mb — > mvM a dipion decay 
Dalitz plot 

Ad = JJdml 2Ma dml^ 2 = jjdsdt (4) 

and pi = p a and Pf = Pb are the three- momenta of the 
initial and final pions (or heavy mesons) in the reac- 
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FIG. 3: ov+t-w+t' (lower) and <VfT'-nr+T ( u PP er ) esti- 
mated from r T /^^.+ 7r - T assuming constant amplitudes. 

tion TiM a — > n Mb in the cm. frame. The Dalitz plot 
for the decay ip' — ► tt + it~ J/ip and the related reactions 
7r + J/ip — > TT + ip' and ir + ip' — ► ir + J/ip are shown as ex- 
amples in Fig.l. The Dalitz plot areas {4o} and dipion 
widths used in this work are given in Table 1; the masses 
assumed are m I+ = 0.1396 GeV, M J/4 , = 3.097 GeV, 
AV = 3 - 686 GeV > m t = 9.460 GeV, M r > = 10.023 GeV 
and My = 10.355 GeV. 

The cross sections for the reactions ir + J/ip — > ir + ip' 
and 7r + '0' — > 7r + J/-0 in the constant amplitude approxi- 
mation are shown in Fig. 2. Evidently the scales of these 
cross sections a few hundred MeV above threshold are ca. 
20 /itb for the endothermic process ir + J/ip — > ir + ip' and 
ca. 0.1 mb for its crossed exothermic partner Tr + ip' — > 
7T+J/-0. This method applied to pion cross sections in 
the upsilon family yields cross sections of ca. 2 /ib for 
7r + T — > 7r + T' and ca. 10 pb for 7r + T' — > 7r + T in the 
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transition 


A c [GeV 4 ] 


IW [keV] 


c [mb GoV 2 ] 


if) ' — » TT + 7V~ J /if) 


0.436 


91.5 ±9.0 


0.65 ±0.06 


T' ^TT + TT-T 


1.023 


8.3 ± 1.3 


0.50 ±0.08 


T" -> TT + TT^T 


6.679 


1.2 ±0.2 


0.012 ±0.001 


T" -> TT + TT^T' 


0.027 


0.7 ±0.2 


1.9 ±0.5 



TABLE I: Experimental dipion-decay Dalitz plot areas, tran- 
sition rates and cross section coefficients. 



analogous kinematic regime (Fig. 3). 



IMPROVED ESTIMATES 

Although the constant amplitude results of the pre- 
vious section are of interest as order-of-magnitude es- 
timates, it is known experimentally that the invariant 
amplitudes for these dipion decays show a strong and 
rather complicated dependence on the tttt invariant mass. 
In ip' — > n + Tr~J/ip a strong suppression of the tttt sys- 
tem at low invariant mass is evident (see for example our 
Fig. 4, or Fig. 6 of Bai et al.\v^). One might have instead 
expected a near-threshold enhancement, due to the at- 
tractive FSI of an S-wave 1=0 7T7t system. In contrast, 
the decay T" — > ttttY has a complicated double-bump 
structure, which does show an enhancement at low tttt 
invariant mass, as well as a high invariant mass enhance- 
ment as seen in ip ' — ► irirJ/ip. This difference between 
the observed tttt distributions suggests that the tttt pro- 
duction amplitude depends significantly on the heavy- 
quarkonium source, and is not accurately described by a 
universal 7T7t amplitude alone. 

Various theoretical models of these dipion decay am- 
plitudes have appeared in the literature. They may be di- 
vided into three broad categories, according to the mech- 
anism that is assumed to give rise to the strong m TT - 
dependence; 1) gluon radiation models 0,112, 2) 
scalar anomaly model s pil I25I Eif . and 3) scalar meson 
exchange models jU l28l l29f. The first two categories 
predict rather similar m T7r dependences, so we will con- 
sider them together. There is also a suggestion that chiral 
symmetry, combined with certain simplifying asumptions 
for amplitudes, can explain the observed m W7r dependence 
in ip' ~ * ttttJ/iP [3(J ; this model suggests an m n „ de- 
pendence similar to the first two categories. Finally, we 
note that several references have considered the decay 
T" — > 7T7rT as a special case, since it has a "double- 
humped" dipion distribution that is not seen in other 
decays 

Gluon radiation and scalar anomaly models assume 
that the important low-energy to T7F dependence is deter- 
mined by aspects of a purely gluonic intermediate state. 
In gluon radiation models this strong energy dependence 
arises from a multipole expansion of the gluon emission 



amplitude, whereas in the scalar anomaly models it is the 
momentum dependence encountered in coupling the glu- 
onic state to the 7T7t final state. Clearly it will be difficult 
to distinguish these two possibilities, although a detailed 
comparison of their predictions with experiment appears 
to favor the scalar anomaly model [l^. 

In the scalar anomaly model one can relate the 7T7r pro- 
duction amplitude to the matrix element of the simplest 
scalar gluon operator between the vacuum and a tttt state, 
(tttt\G^G^\0) . This matrix element can be determined 
because the operator G a G a is proportional to the tri- 
angle anomaly in the trace of the energy-momentum ten- 
sor, which in a low energy pion effective lagrangian is 
quadratic in the pion field. The matrix element of this 
operator gives a near-threshold dependence of A cx m%„, 
which is A cx t in our kinematics. 

Meson exchange models assume that the tttt system is 
produced by an intermediate scalar /o meson (often re- 
ferred to as the "er" ) , and the observed m^-dependence 
at higher invariant mass is due primarily to this meson. 
Fits to the dipion decay data (excluding the problemati- 
cal T" — > 7T7rT) typically prefer a light, broad state with 
a mass near 0.5 GeV; see for example [28l 123 . l33 j . 

Ideally we would cross an accurate model of the 7T7t 
production amplitudes into the tt scattering regime to 
estimate tt + (QQ) closed-flavor cross sections, but this 
is not yet possible because no model gives a good simul- 
taneous description of all the experimentally observed cc 
and bb dipion mass distributions. For the present we 
will simply assume the near-threshold dependences 
suggested by the existing models, and evaluate the cross 



sections predicted for tt + J/rp 



tt 



1b 1 in each case. 



For the gluon radiation and scalar anomaly models we 
assume 



(\A(s,t)\ 2 ) =c 2 t 2 



(•5) 



For a power-law form (|.4(s,i)| 2 ) = c n t n ) the relation 
between decays and cross sections Eq.(3) generalizes as 
follows; the decay rate Eq.(2) becomes 
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where I^ n ' is the integral of t n over the A4b — > TnrM. a 
Dalitz plot, 



ds dt t n 
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The constant c n is determined by the measured dipion 
partial width using Eq. (6) , and substitution of c n t n into 
the cross section formula Eq.(l) then gives 
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K(\h\ n+1 -\h\ n+1 ) . (8) 



We again specialize to the process ir + J/ip — » 7r + V'- 
Setting n — 2 (scalar anomaly models), for Eq.(7) we 
find Z( 2 ) = 0.01954 GeV 8 . The cross section in Eq.(8) 
may then be evaluated numerically, which gives the result 
shown in Fig. 5. Note that the scalar anomaly model leads 
to a rapid increase of the cross section relative to the 
constant amplitude model above s/s « 3.9 GeV. This is 
due to the i 2 -weighting combined with the rapid increase 
in the range of t covered by this reaction with increasing 
y/s, which is evident in Fig.l (lower right region). 

Finally, for our meson exchange model we assume a 
generalized Breit-Wigner form which incorporates the 
scalar anomaly soft-pion factor. For t > this is 



(\A(s,t)\ 2 ) = 
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C 2 t 



(Vi-M fo r + T%/4 



(9) 



where Mf and T / are the mass and width of the hypo- 
thetical scalar meson source of the irir events. We have 
incorporated the t 2 scalar anomaly soft-pion dependence 
in Eq.(9) because pure Breit-Wigner forms required an 
unrealistically narrow f (T f w 100-150 MeV) and gave 
rather poor fits to the data. In contrast the hybrid form 
Eq.(9) clearly gives an acceptable fit (Fig. 4), although 
we emphasize that this meson exchange model is unphys- 
ical because the fitted parameters Mf — 536 MeV and 
T f — 260 MeV are inconsistent with the experimental 
1=0 7T7T S-wave phase shift. 

We can again use Eqs.(l,2) to determine the cross sec- 
tion for 7r + J/-0 — > ir + ijj' implied by this decay model 
(ignoring the problem of disagreement with phase shifts) . 
The result is 
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FIG. 4: Fitting a meson exchange model to the tt-k distribu- 
tion observed in /_,„.+„.- j/^ . 
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(10) 

Note that this integration is over negative values of t, 
whereas the decay rate integral is over positive t, which 
leads to different signs in the Breit-Wigner functions. 
The Dalitz plot decay integral 



2>„ = 



ds dt 



t 1 



(Vt-M fo r+Tl/4 



(11) 



for ip' — ► ir+Tr^J/tp equals 0.08069 GeV 6 given our ex- 
ternal meson masses and fitted /o parameters. The inte- 
gral over t in Eq.(10) is 
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FIG. 5: The near-threshold cross section for -k + J/ip — > -K + ip ', 

in constant amplitude, scalar 



estimated from ^^'-, w + n - j/^, 
anomaly and meson exchange models. 



(2 - 20c 2 + 10c 4 ) 



=i+V=tT/M /0 



(12) 



where c = Tf /2Mf . Combining Eqs. (10-12) gives 
the meson exchange model prediction for j/^^^+^p i , 
which is also shown in Fig. 5. Evidently the meson ex- 
change cross section is suppressed near threshold, due to 
the separation from the fo pole. 



COMPARISON WITH PREVIOUS WORK 

Two of the theoretical references on dipion decays cited 
previouslv, Sorge, Shuryak and Zahed [22] and Chen and 
Savage [23, also discussed results for closed- flavor pion 
scattering cross sections with heavy quarkonia, specifi- 



cally for n + J/tp — ► n + ip' at low energies. Since these 
references actually used dipion decay data to normalize 
their scattering amplitudes, approximate agreement with 
our results should be anticipated. Fujii and Kharzeev [34| 
have also evaluated closed-flavor charmonium cross sec- 
tions, using a color-dipole scattering model. They do not 
use dipion decay data as a direct input, although they do 
note that their two results for the rate T(ip' — > nnJ/ip) 
are not far from experiment. 

In the earliest reference , Sorge et at assume scalar 
meson exchange with a high-mass /o(1400) as the nn 
source. (They make the important and often neglected 
observation that assuming a low-mass /o of only moder- 
ate width, as in Fig. 4, disagrees with the experimental 
1=0 nn S-wave phase shift.) Their cross section close to 
threshold is shown in Fig. 6, and is evidently qualitatively 
similar to our scalar anomaly model result. Chen and 
Savage |2^| used a gluon radiation model to describe this 
reaction, and quoted cross sections at tree level and with 
one-loop chiral corrections. These results are also shown 
in Fig. 6, and are numerically rather similar to Sorge et al. 
Finally, Fujii and Kharzeev [34| used a color-dipole scat- 
tering model, and quote results for this process both with 
and without a scalar /o form factor for the nn system. 
(Their form factor is inferred from the 1=0 nn S-wave 
phase shift.) Their results without a form factor are sim- 
ilar to the earlier predictions. With a form factor they 
find a much smaller cross section, which is rather close 
to our scalar meson exchange result. 

SUMMARY AND CONCLUSIONS 

The search for the quark gluon plasma has led to great 
interest in the scale of the cross sections of heavy quarko- 
nia interacting with light hadrons near threshold. Unfor- 
tunately, little is known about these cross sections ex- 
perimentally. In this paper we have used crossing sym- 
metry and several simple amplitude models to estimate 
the closed-flavor cross sections for heavy quarkonia scat- 
tering on pions near threshold, using the experimentally 
known dipion decays as input. The method is applied 
both to charmonia and to the bb system. For the sim- 
plest cases of IS <-> 25 transitions, assuming constant 
amplitudes we estimate the cross sections a few hundred 
MeV above threshold to be ca. 20 fxh for nJ/ip — ► nip', 
and ca. 2 /ib for nT — > nT'. The corresponding time- 
reversed, exothermic reactions are estimated to be about 
0.1 mb and 10 /xb respectively. 

We note that the strong dependence of the decay am- 
plitude out = observed experimentally in the dipion 
decays near threshold makes analytic continuation to the 
pion scattering regime rather problematic. The depen- 
dence of the cross section for n J/ ip — >• nip ' (used as our 
example) on the different model amplitudes is clearly ev- 
ident even near threshold (Fig. 5). For this reaction the 
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FIG. 6: A comparison of our theoretical cross sections for 
n + J/%p — > -K + ip' (from Fig. 5) with the results of Sorge, 
Shuryak and Zahed [27) . Chen and Savage [23f and Fujii and 
Kharzeev [3^1 . Note that of the three forms we have assumed, 
only the lowest cross section (meson exchange model) gives a 
good fit to the observed nn distribution (Fig.4). 

various models we considered gave consistent cross sec- 
tions of ~ 5-15 /ib at E cm = 3.9 GeV, but the predictions 
diverged rapidly with increasing invariant mass. 

In future, high statistics studies of dipion decays at 
CLEO might provide additional useful information about 
the decay amplitude. In particular, it would be useful 
to accurately determine the s dependence of the decay 
amplitude in the Dalitz plot of Fig.l experimentally, in 
addition to the already well known t = rri^ n dependence. 

We note in passing that the cross sections for these 
near-threshold closed-flavor processes are much smaller 
than the millibarn scale typically found for open-flavor 
reactions such as n J/ip -> D*D 0, 0, H H 011 ©, OH- 
The possibility that the rather weak closed-flavor reac- 
tions are due to open-flavor scattering at second order 
|3lL I35L l36| , rather than to exchange of an unphysically 
light scalar meson, is an interesting suggestion which 
merits future investigation. 
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